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The masses of two-gluon glueballs are studied with a semirelativistic potential model whose inter- 
action is a scalar linear confinement supplemented by a one-gluon exchange mechanism. The gluon 
is massless but the leading corrections of the dominant part of the Hamiltonian are expressed in 
terms of a state dependent constituent gluon mass. The Hamiltonian depends only on 3 parameters: 
the strong coupling constant, the string tension, and a gluon size which removes all singularities 
in the leading corrections of the potential. Accurate numerical calculations are performed with a 
Lagrange mesh method. The masses predicted are in rather good agreement with lattice results and 
with some experimental glueball candidates. 
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I. INTRODUCTION 



o 
o 

(N 

Q I The existence of bound states of gluons, called glueballs, is a prediction of the QCD theory. The experimental 
D ■ discovery of such particles would give a supplementary strong support to this theory. But, a reliable experimental 
identification of glueballs is difficult to obtain, mainly because glueball states might possibly mix strongly with nearby 
meson states. Nevertheless, the computation of pure gluon glueballs remains an interesting task. This could guide 
experimental searches and provide some calibration for more realistic models of glueballs. 

The potential model, which is so successful to describe bound states of quarks, is also a possible approach to study 
y—i glueballs. Among the pioneer works using this formalism, the one of Corwall and Soni is particularly interesting 
^ ' Assuming a nonrelativistic kinematics and a saturated confinement supplemented by a one-gluon exchange interaction, 
masses of pure gluon glueballs were computed. However, only the four lightest two-gluon states (0++, +, 1 2++) 
[ were computed and found between 1.2 and 1.8 GeV. Using a similar model, two-gluon glueballs have been studied in 
. Ref. 0. The masses of states with L = 0, 1 and 2 orbital momentum were computed and several states are found 
I ■ below 3 GeV. Unfortunately, we think that this model suffers from several drawbacks which spoil any possible physical 
, conclusions. This has been discussed elsewhere 0- Nevertheless, we think that these drawbacks can be corrected in 
' order to obtain a more reliable potential model. 

In this paper, we compute two-gluon glueball masses using various modifications of the potential model obtained 
^ ^ two decades ago by Cornwall and Soni. After a critical study of these various models, we conclude that a spectra in 
I rather good agreement with lattice calculations and some experimental glueball candidates can be obtained, provided 
several conditions are fulfill: a semirelativistic Hamiltonian is used, the gluon has a finite size, the confinement is a 
scalar interaction, and a dynamical constituent gluon mass is used in the leading relativistic corrections. 

In Sec. m the Hamiltonian model is presented with the notion of gluon size. Three variants of our model are 
discussed in Sec. Illll and a glueball spectrum is presented. Concluding remarks are given in Sec. II VI and some useful 
technical details are given in appendix. 



II. HAMILTONIAN 



The two-gluon Hamiltonian contains a kinetic part Hq, a short-range part Vsr due to the one-gluon exchange 
between the two valence gluons, and a confining interaction V^onf, as the model proposed by Corwall and Soni Q 

H^Ho + VsR + Konf. (1) 

Following Refs. there is no constant potential, contrary to usual Hamiltonians for mesons and baryons. This 

model can be considered with both nonrelativistic (Schrodinger equation) and semirelativistic (spinless Salpeter 
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equation) kinematics [fi ^ c ^ 1) 



Ho = 2mK + — or Hq = 2Wp2 + rn\, (2) 

TTlK ^ 

where rriK is the effective gluon mass appearing in the free part of the Hamifronian. If a Schrodinger Hamiltonian is 
used, it is necessary to verify that, for each glueball state, the quantity {p"^ /mj^), which can be considered as the 
mean speed of a gluon, is small with respect to 1. 



A. Short range potential 

We use the short-range potential between two gluons proposed in Ref. I*]. After some manipulations, this potential 
takes de following form 



2s - Im^ 1 
Qw? 3 
Att ,o , , / - 2s 5 



nr) ^^^^^ + ^t;S' 
L S 



3A U'{r) 



2m? 

+ A ^ith U{r)^- , (3) 

where L and S are the usual orbital momentum and spin operators, and where 

T ^ S'^ ^ iiS ■ rf (4) 

is the tensor operator, m is an effective gluon mass, which can differ from the mass itlk (see below). The quantity 
s is in principle the square of the glueball mass, but we will always take s — Am? as it is suggested in Ref. X\. The 
parameter A is linked to the strong coupling constant as by the relation |5(| 

A = 3as. (5) 

This potential has a priori a very serious flaw: depending on the spin state, the short distance singular parts of the 
potential may be attractive and lead to a Hamiltonian unbounded from below. We will see in Sec. Ill CI how to cure 
this problem. 



B. Confinement potential 

The dominant part of the interaction between the two gluons is the confinement. As the leading relativistic 
corrections are taken into account in the short-range part of the interaction, it is natural to keep the same order 
corrections for the confinement potential. The Lorentz structure of the confining interaction is not well known yet. 
In this work, we follow the prescription of Ref. 0: if the radial form (static or zero order part) of the confinement is 
Wx {f) , then the total confinement interaction is written 

Fconf = Wxir) - -l^W^ir) L ■ S, (6) 
2m^r 

where the effective mass m is the same than the one appearing in potential Q • Actually, this form contains only the 
dominant correction, which is a spin-orbit contribution, and it is only valid for large values of the distance r. But, 
for our purpose, this approximation is sufficient. The interaction ((HJ corresponds to a confinement with a dominant 
scalar structure Let us note that the spin-orbit contribution from the confinement counteracts the spin-orbit 
contribution from the one-gluon exchange, and plays an important role to obtain a spectra in agreement with lattice 
calculations. 

Two radial forms Wx{r) can be used. In Refs. 0,0,13, the confinement potential saturates at large distances 



Wp{r) = 2m (1 - e-'^"'') 



(7) 



Such a form can simulate the breaking of the color flux tube between the two gluons due to color screening effects. 
The maximal mass for a glueball is then 4m. Another simpler form is proposed in Ref. 

Wa{r)^aGr = ^ar, (8) 

where ac is the string tension between two gluons and a the usual string tension between a quark and an antiquark. 
The 9/4 factor is due to the color configuration of the gluons These two potentials coincide at small distances, 

which implies that 

Potential ((SJ seems a priori inappropriate since strings joining gluons must always break if a sufficiently high energy 
is reached. But this phenomenon must only contribute to the masses of the highest glueball states. 



C. Gluon size 



Within the framework of a potential QCD model, it is natural to assume that a gluon is not a pure pointlike particle 
but an effective degree of freedom which is dressed by a gluon and quark-antiquark pair cloud. Such hypothesis for 
quarks leads to very good results in the meson and baryon sectors. We assume here a Yukawa color charge 
density for the gluon 

= ^ ' (10) 

where 7 is the gluon size parameter. The interaction between two gluons is then modified by this density, a bare 
potential V being transformed into a dressed potential V. This potential is obtained by a double convolution over 
the densities of each interacting gluon and the potential, ft can be shown that this procedure is equivalent to the 
following calculation 

V{r)^ I dr'V{r')T{r-r') with r(M) = -^e""/'^. (If) 
J 87r7'^ 

Convolutions for some useful potentials are given in appendix 1X1 

Other color charge densities could be used, a Gaussian one for instance We have nevertheless strong indications 
that such a change cannot noticeably modify the results jf^l . We choose the Yukawa density because all convolutions 
are analytical with this form. 

The convolution (|f f |) of potential e~™'^ /r with the color density (|f 0|) removes all singularities in the short-range 
interaction For consistency, we apply the same regularization to the confinement potential ((H)), although no 
singularity are present in the radial form Wx{r). 



III. RESULTS 



A. Some general considerations 



Nonrelativistic potential models have been intensively used to compute static properties of mesons and bar yon s. 
But numerous works show that semirelativistic potential models can give better results (see for instance Refs. |f 3lll4| ). 
Although the gluon effective mass is expected around 700 MeV |^ , which is heavier than the assumed constituent 
strange quark mass, the relevance of a nonrelativistic dynamics for the gluon is questionable. Using the various models 
discussed below with several different sets of realistic parameters, we have always obtained values of V\PV™1T around 
unity (and sometimes largely above) when a nonrelativistic kinematics was used (if mn ^0). As we shall see below, 
the model III is by nature a semirelativistic one. Both kinematics can be used for models I and 11, but we have verified 
that the drawbacks of these two models cannot be solved by a change of kinematics. Consequently, in the following, 
we will only present results from spinless Salpeter Hamiltonians. 

In order to avoid singularities, potentials U{r) and Wx{r) of formulas © and (0 have been replaced by their 
corresponding dressed forms U{r) and Wx{r). Another way to get rid of singularities in the short-range potential 
is to treat Vsr as a perturbation (at least when it is attractive) of the dominant confinement interaction. But, for 



a lot of states, computed with different sets of realistic parameters, the contribution of the short-range part can be 
comparable to the one of the confinement part. So, a perturbative approach of the singularities can hardly be justified. 
In the following, we will always treat the potential Vsb. non perturbatively. 

The tensor operator T is responsible of channel coupling. Its matrix elements are given in appendix It can be 
shown that the total spin S of two gluons is always a good quantum numbers, but mixing of orbital momenta with 
the same parity is possible. For instance, the 2++ glueball with 5 = 2 is the mixing of three states with L — 0, 2 
and 4. But, it is not coupled with the 2++ glueball with 5 — and L — 2. In principle, the effect of mixing is a 
second order correction with respect to the contribution of the diagonal term. In this paper, results are only shown 
when off-diagonal tensor contributions are neglected. Nevertheless, we will discuss below the effect of mixing for our 
various models. 

The general characteristics of all our models are: 

• Semirelativistic kinematics; 

• All radial forms convoluted following relation Hll|l: 

• VsR not considered as a perturbation of Vconf] 

• No channel coupling with T. 

The eigenvalue problem has been solved by the Lagrange-mesh method which allows a great accuracy as well as for 
Schrodinger equation [15] as for spinless Salpeter equation ll^. 

In the models considered below, the gluon masses {mn and m) will be fixed by physical considerations. We are 
then left with 3 parameters: as (or A) for Vsr, (3 or a for Vconf, and the gluon size 7 for which less constraints exist on 
its value. Fortunately, the mass of the lightest 2++ state is nearly independent of the parameters ag and 7. For this 
state, the spin-orbit and diagonal tensor potentials vanish and the two remaining contributions have opposite signs. 
So, the confinement potential is the largely dominant contribution to this 2++ mass (see Fiars. 1113)1 . We can fix the 
value of /3 or a with this state only, knowing that a lattice calculation and a quasiparticle model favor a 2+"'" 
mass around 2.4 GeV, but that some experimental candidates are found around 2 GeV 0,1131 • Some agreements 
between theoretical calculations and experimental data exist about mass ratios of some lightest glueball candidates 
(see Refs. [Ulli |20| and TablelD: M{0++)/M{2++) - 0.72 - 0.78 and M{0-+)/M{2++) ~ 1.10. So, we wiU fix 
the parameters as and 7 in order to reproduce at best these mass ratios. Let us note that the mass ratios in lattice 
calculations can be more interesting quantities to consider than absolute masses, due to the existence of normalization 
problems jl7i |. 



B. Model I 

The first model we consider is by two aspects close to the models of Refs. Q and 0: The short-range part is 
supplemented by a saturated confinement potential and the gluon is assumed to be characterized by an unique effective 
mass rriK = rn, around the typical value of 700 MeV. In this paper, we choose the value of Ref. Nevertheless, 
the model is semirelativistic and the short-range part is not treated perturbatively. The particular characteristics of 
model I are: 

• rriK = m = 0.670 GeV 0; 

• Konf with Wi3{r). 

To find a mass of the lightest 2++ state around 2 GeV, it is necessary to take f3 « 0.2, which corresponds to 
a w 0.08 GeV^, a quite unrealistic value. A mass around 2.4 GeV can be obtained with f3 « 0.5 corresponding to 
a « 0.2 GeV^, a more realistic value. The parameters (3 — 0.5, as — 0.5, and 7 = 0.5 GeV^^ give the following 
lightest mass ratios: M(0++)/M(2++) = 0.93 and M(0-+)/M(2++) = 0.73. Let us note that if the spin-orbit 
contribution from confinement is not taken into account, the last ratio decreases to 0.61. A lower value for the mass 
ratio M(0"'""'")/M(2"'"+) can be obtained by modifying the parameters as and 7, but in this case, the mass ratio 
M(0 ^)/M(2'^~^) also decreases. It can even take negative values for realistic values of the parameters as and 7. 
This nonphysical behavior is due to the spin-orbit potential from Vsr which becomes very attractive. It is worth 
noting that a variation of the gluon effective mass of 200 MeV around the value 670 MeV does not change noticeably 
the results. Thus, the model I is not able to describe neither the lightest 2++ mass state, nor the lightest mass ratios 
M(0++)/M(2++) and A/(0-+)/M(2++). 

If the channel mixing due the tensor operator is turned on, the situation gets worse. For instance, without channel 
mixing the lightest 2++ state {L = 0, S — 2) has a reasonable mass. With channel mixing, the L — Q component is 
coupled with L = 2 and L ~ A components for which the total spin-orbit contribution is very attractive. The mass of 
this state becomes then negative, even for realistic values of the parameters /?, as and 7. 



C. Model II 



The model II is the same as model I but with the saturated confinement replaced by the linear confinement. The 
particular characteristics of model II are: 

• rriK = m = 0.670 GeV 0; 

• Konf with Wair). 

This model and the previous one give essentially the same results about the lightest mass ratios Af (0++)/Af (2++) 
and M(0 ^)/Af(2++). Moreover, masses around 2.4 GeV and 2 GeV are obtained for the lightest 2++ state with 
a Ki 0.12 GeV^ and a « 0.07 GeV^ respectively, quite unrealistic values of the meson string tension. 

D. Model III 

With the two previous models, the spin-orbit effect from the one gluon-exchange is too attractive and cannot be 
counteracted efficiently by the spin-orbit contribution coming from the confinement. The strength of this attractive 
potential can be reduced by decreasing the values of as- But in this case, it is not possible to obtain reasonable mass 
ratios for the lightest 0++, 2++, and ^ states. Another possibility is to increase the values of the effective mass m. 
But, if we keep the link ttik — m, then too high masses are obtained for all glueballs, due to the contribution of the 
kinetic part of the Hamiltonian. Fortunately, it is physically relevant to choose mn <SC m. 

For a system of two identical particles with mass m, the coefficient 1/rn? appears naturally in the relativistic 
corrections of a static potential. X better approximation, proposed in Ref. and used for instance in Ref. j21j| , is to 
replace this coefficient by 1/E'^{p) where E{p) = ^/p^ + rn^ . 

A similar procedure is also proposed within the auxiliary field formalism (also called einbein field formalism) |22j |. 
which can be considered as an approximate way to handle semirelativistic Hamiltonians 23]. Within this approach, 
the effective QCD Hamiltonian for two identical particles (quark or gluon) depends on the current particle mass m 
and on a state dependent constituent mass /i = (p^ + m^) 23] . All corrections to the static potential are then 
expanded in powers of 

We will adapt these prescriptions in the model III. In principle, the mass appearing in the leading corrections must 
be replaced by the operator E{p) = \/p^ + vn? where m is the mass appearing in the kinetic operator. This leads to 
a very complicated nonlocal potential which is very difficult to handle. So we will use an approximation. 

Following the hypotheses of Ref. , we will assume that the gluons are massless and that the dominant effective 
QCD Hamiltonian Hgg for a two-gluon glueball is written 

H<,g = 2v/p2 + acr. (12) 

The eigenvalues M„i and the constituent masses /i„L corresponding to this Hamiltonian are in very good approxi- 
mation given by the following relations [23l l23 ] 

/g,^^\ 3/4 

M„i, = 4/i„L with y^riL = (~|~ j ■ (13) 

e„L is an eigenvalue of the Hamiltonian -t- \x\, in which q and x are dimensionless conjugate variables. These 
eigenvalues can be accurately computed with a Lagrange- mesh method for instance 15] . Let us note that e„o is the 
nth zero of the Airy function. These constituents masses will then appear in the leading corrections to the dominant 
hamiltonian Hgg. 

Finally, the particular characteristics of model HI are: 

• mK — and ra — ^nL (I13|l : 

• V'conf with Wa{r). 

Let us note that, with our hypotheses, the constituent mass depends on the principal quantum number n of the 
glueball, on its orbital momentum L, but not on the quantum numbers S and J. 

As mentioned before, we fix the value of the string tension only with the mass of the hghtest 2++ glueball. 
Equations (|13|l show that, in first approximation, the mass scale is simply given by ^Joq. By computing a great number 
of spectra for various parameters, we have remarked that the values of the lightest mass ratios M (0++)/Af(2++) and 
M(0 ^)/Af(2++) cannot be fixed independently. Provided an approximate linear dependence is kept between the 



two parameters a and 7, these mass ratios do not change significantly. Finally, we have chosen to present the glueball 
spectra for two sets of parameters. 

With a = 0.16 GeV^ as = 0.40, and 7 = 0.504 GeV^^ (set A), the mass of the lightest 2++ gluebaU is 2051 MeV, 
which is close to some experimental candidates. These values for a and as are near those used in some recent baryon 
calculations I2H. Moreover, the value of the string tension is close to a value found in a recent lattice study [26l |. 
With a = 0.21 GeV^, as = 0.50, and 7 = 0.495 GeV~i (set B), the mass of the lightest 2++ glueball is 2384 MeV, 
which is close to a result obtained with lattice calculations. All results are presented in table HI and in Fig. 01 with 
some results from a lattice calculations ^| and from a quasiparticle approach with no free parameters jlgj , and with 
some possible experimental candidates |l9ll20l|. 

We can see that the mass ratios for the two sets of parameters are in rather good agreement with the theoretical 
mass ratios predicted by the lattice study Moreover, the absolute masses for set B are within the theoretical 

error bars of the lattice masses. The largest discrepancy is for the first excited 0++ glueball, the state predicted by 
the lattice calculations with the largest error bar. 

Our mass ratios are similar to those obtained in the quasiparticle model (ISj . but are closer to the mass ratios of 
the lattice studies. It is worth mentioning that this quasiparticle approach contains no free parameters. Again, it 
favors a 2.4 GeV value for the mass of the lightest 2++ state, like the lattice model. 

The lack of reliable identification of glueballs makes comparison with experimental data more hazardous. The set 
A results are in rather good agreement with data only for the lightest 0++, 2++, and glueball candidates, the 
states used to fix the parameters. The general tendance of our model and of lattice calculations is to predict excited 
states with higher masses than those which seem to be observed. 

Lattice calculations seem to rule out the presence of 1 ^ and 1++ states below 4 GeV. This can be qualitatively 
understood in terms of interpolating operators of minimal dimension which can create glueball states, with the 
expectation that higher dimensional operators create higher mass states: the lowest states 0++, 2++, ^ and 2 ^ 
are produced by dimension-4 operators, while 1^"*" and 1 are respectively produced by dimension-5 and dimension-6 
operators JJ}. Nevertheless, our model predicts the existence of 1 ^ and 1++ states around 3 GeV. In the model of 
Ref. 0, a 1 ^ state is predicted below the 2++ state (no 1++ state is mentioned). Possible experimental candidates 
exist for low mass 1++ states but as mentioned above, the identification is far from certain. The presence of 
the relatively low mass 1 and l"*""*" states in our model may be due to the use of massive valence gluons with three 
states of polarization (creation of a spin one glueball with two massless gluons, with only two states of polarization, is 
problematical). In our model III, the gluon is massless in the kinetic part, but a constituent non-zero mass unavoidably 
appears for the spin corrections Q . The presence of spin one states around 3 GeV in our model could indicate what 
are the limits of a potential approach. 

If the spin-orbit contribution from confinement is not taken into account, the agreement between our masses and 
the lattice results become poorer. For the parameters of set A, the mass ratio for the lightest ^ glueball changes 
from 1.06 to 0.87, and the mass ratio for the first excited 0++ glueball changes from 1.26 to 0.83. This shows that 
the spin-orbit contribution from confinement is an important ingredient of the model. 

The channel mixing due to the tensor force is difficult to implement within this model. As fi depends on the orbital 
momentum, the diagonal potential for each channel is characterized by a different value of /i. The problem is to define 
this parameter for the mixing potentials. We have performed several test computations using a mean value of fi for all 
channels. This gave us strong indications that the coupling of channels has a small influence on the glueball masses, 
contrary to the two previous models. We estimate that the masses of the lightest glueballs could be modified by a 
quantity comprised between 50 and 100 MeV. 



IV. CONCLUSION 



The masses of pure two-gluon glueballs have been studied with a semirelativistic potential model. The potential 
is the sum of a one-gluon exchange interaction and a linear confining potential, assumed to be of scalar type. The 
gluon is massless but the leading corrections of the dominant part of the Hamiltonian are expressed in terms of a 
state dependent constituent mass. The Hamiltonian depends only on 3 parameters: the strong coupling constant, 
the string tension, and the gluon size. This last parameter, less constrained than the two others by the QCD theory, 
removes all singularities in the leading corrections of the potential. These corrections are not treated as perturbations 
of the dominant part. All masses have been accurately computed with a Lagrange mesh method. 

The masses predicted by our potential model are in agreement with experimental glueball candidates only for the 
lightest 0++ 2++, and 0"+ states [Tl l20j. but are in rather good agreement with spectra obtained by a lattice 
calculation and in reasonable agreement with spectra obtained by a quasiparticle model . A notable difference 
is the presence in our model of spin one states around 3 GeV. This could indicate the limit of the validity for a 
potential approach. 



We have tested other nonrelativistic and semirelativistic potential models in which a constant constituent gluon 
mass is used, and we have found that it is not possible to obtain good spectra for realistic values of the QCD parameters 
(see Sec. lmBl and lmCl . The main problem arises from the strongly attractive spin-orbit potential for the one-gluon 
exchange. When its strength is not reduced by a large constituent gluon mass, it can lead to negative nonphysical 
glueball mass. 

The constituent gluon mass is introduced in our model by an approximate procedure which relies on the existence 
of a pure linear confinement between the gluons [2^1 123 | . A more physical ansatz should be to define the constituent 
mass as a momentum dependent operator for massless gluon) 7]. It could then be possible to take into account 

correctly the channel coupling due to the tensor forces, and to use naturally a saturated confinement potential. It 
could also be interesting to compute three-gluon glueball masses within the same model. Such a work is in progress. 
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APPENDIX A: CONVOLUTIONS 

Applied for some useful potentials, the formula Hll|l gives 

P{r) = ^e-^/\ (Al) 

F = r + ^ f 1 - e-''/A - 76-''/^, (A2) 



3 — ar 



r 



(0^7^ — 1)3 \ ^ J (0^7^ — 1)2 

1 /e-'"' e-''/^\ e-''/^ 



(A3) 

One can easily verify that lim^^o U{r) — U{r) for each potential U{r). 

APPENDIX B: ANGULAR MOMENTUM OPERATORS 

A system of two particles, with spin si and S2 respectively, with a total spin S and a total orbital angular momentum 
L coupled to a total angular momentum J, is noted here \L^S) = \si, S2] S] J) . The mean value of the operators 
52 and L ■ S are trivial to compute 

{L\S'\S''\L,S) = S{S + 1)6l'^l6s',s. (Bl) 

{L',S'\L- S\L,S) = ]^[j{J+1)-L{L+1)~S{S+1)\5l.l'5s,s'- (B2) 

The computation of the mean value of the operator T is much more involved. Using formulas from Ref. |27| , one can 
find {fi = ^/2n + 1) 

{L'^S'\T\L^S) = {^if+^'+^'+^SS'LV f i n ( I (B3) 
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our model] obtained with model HI for two sets of parameters (A: a = 0.16 GeV^, as = 0.40, and 7 — 0.504 GeV~^; B: 
a = 0.21 GeV^, as ~ 0.50, and 7 = 0.495 GeV^^). Some theoretical results from other models and some possible experimental 
candidates are also indicated. The lightest 0^^, 2"*"+, and 0""*" states are taken as inputs to fix the parameters. 





[L,S\ 


Model HI 
A 


B 


Lattice [171 


Quasiparticle 


Experiment 


Experiment [2^ 


0++ 


[0,0] 


1604 


(0.78) 


1855 


(0.78) 


1730±50±80 (0.72) 


1980 


(0.82) 




1507±5 (0.78) 




[2,2] 


2592 


(1.26) 


2992 


(1.26) 


2670±180±130 (1.11) 


3260 


(1.35) 




2105±15 (1.09) 




[0,0] 


2814 


(1.37) 


3251 


(1.36) 












2++ 


[0,2] 


2051 


(1.00) 


2384 


(1.00) 


2400±25±120 (1.00) 


2420 


(1.00) 


2020±50 (1.00) 


1934±12 (1.00) 




[0,2] 


2985 


(1.46) 


3447 


(1.45) 




3110 


(1.29) 


2240±40 (1.11) 






[2,0] 


3131 


(1.53) 


3611 


(1.51) 








2370±50 (1.17) 






[2,2] 


3230 


(1.57) 


3695 


(1.55) 












0-+ 


[1,1] 


2172 


(1.06) 


2492 


(1.05) 


2590±40±130 (1.08) 


2220 


(0.92) 


2140±30 (1.06) 


2190±50 (1.13) 




[1,1] 


3228 


(1.57) 


3714 


(1.56) 


3640±60±180 (1.52) 


3430 


(1.42) 






1-+ 


[1,1] 


2626 


(1.28) 


3011 


(1.26) 














[1,1] 


3349 


(1.63) 


3852 


(1.62) 












2-+ 


[1,1] 


2573 


(1.25) 


2984 


(1.25) 


3100±30±150 (1.29) 


3090 


(1.28) 


2040±40 (1.01) 






[1,1] 


3345 


(1.63) 


3862 


(1.62) 


3890±40±190 (1.62) 


4130 


(1.71) 


2300±40 (1.14) 




1++ 


[2,2] 


3098 


(1.51) 


3501 


(1.47) 








~ 1700 (0.84) 






[2,2] 


3753 


(1.83) 


4294 


(1.80) 








2340±40 (1.16) 




3++ 


[2,2] 


3132 


(1.53) 


3611 


(1.51) 


3690±40±180 (1.54) 


3330 


(1.38) 


2000±40 (0.99) 






[2,2] 


3762 


(1.83) 


4332 


(1.82) 




4290 


(1.77) 


2280±30 (1.13) 




4++ 


[2,2] 


2897 


(1.41) 


3360 


(1.41) 




3990 


(1.65) 


2044±? (1.01) 






[2,2] 


3633 


(1.77) 


4197 


(1.76) 




4280 


(1.77) 


2320±30 (1.15) 
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FIG. 1: Masses of glueball states 0++ {L = S = 0), 2++ {L = 0,S = 2), and 0"+ (L = S = 1) as a function of the string 
tension a, for the model III with as = 0.5 and 7 = 0.52 GeV~^. 
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FIG. 2: Masses of glueball states 0++ (L = 5 = 0), 2++ (L = 0, S = 2), and 0"+ (L = 5 = 1) as a function of the gluon size 
7, for the model III with as = 0.5 and a = 0.2 GeV^. 
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FIG. 3: Masses of glueball states 0++ (L = 5" = 0), 2++ (L = 0, 5 = 2), and 0~+ = 5 = 1) as a function of the strong 
coupling constant as, for the model III with a = 0.2 GeV^ and 7 — 0.52 GeV~^. 
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FIG. 4: Glueball mass ratios normalized to the lightest 2"*"^ sta te ( see Table^l. Cross: Set of parameters B for model III; Black 
circle: Lattice results fl^ : Black square: Quasiparticle model llSlI : White circle: Experiment p^: White square: Experiment 
[2nl |. The experimental states with a question mark are seen but the uncertainty is not known. 



